On Computational Power of Quantum 
Read-Once Branching Programs 



Farid Ablayev* 

Kazan Federal University, 
Institute for Informatics of Tatarstan Academy of Sciences 
Kazan, Russian Federation 

f ablayevOgmail . com 

Alexander Vasiliev ^ 

Kazan Federal University, 
Institute for Informatics of Tatarstan Academy of Sciences 
Kazan, Russian Federation 

Alexander . KSUOgmail . com 

In this paper we review our current results concerning the computational power of quantum read-once 
branching programs. First of all, based on the circuit presentation of quantum branching programs 
and our variant of quantum fingerprinting technique, we show that any Boolean function with lin- 
ear polynomial presentation can be computed by a quantum read-once branching program using a 
relatively small (usually logarithmic in the size of input) number of qubits. Then we show that the 
described class of Boolean functions is closed under the polynomial projections. 

I Introduction 

One notable thing about the recent realizations of a quantum computer (say, the one based on multiatomic 
ensembles in resonator ifTTI ) is that they abide the isolation of a quantum system, e.g. in [11] the trans- 
formation of a quantum state is performed by an external magnetic field. Thus it is quite adequate to 
describe such computations by quantum models with classical control. Here we consider one of such 
models - the model of quantum branching programs introduced by Ablayev, Gainutdinova, Karpinski 

I I ] {leveled programs), and by Nakanishi, Hamaguchi, Kashiwabara fl4l {non-leveled programs). Later 
it was shown by Sauerhoff |[T6l that these two models are polynomially equivalent. 

It is also worth noting that in spite of constant progress in experimental quantum computation all 
of the physical implementations of a quantum computer are still rather weak in a sense that they are 
suffering from fast decoherence of the quantum states and are able to organize the interaction of a small 
number of qubits. This naturally leads to the restricted variants of a quantum computer - the idea first 
proposed by Ambainis and Freivalds in 1998 [6]. Considering one-way quantum finite automata, they 
suggested that the first quantum-mechanical computers would consist of a comparatively simple and fast 
quantum-mechanical part connected to a classical computer. 

In this paper we consider a restricted model of computation known as Quantum Read-Once Branch- 
ing Programs of polynomial width. The classical variant of this model is also known in computer science 
as Ordered Binary Decision Diagrams (OBDDs) and that is why we will also use the notion of quantum 
OBDDs (QOBDDs) for the considered model. The small coherence time is formalized in this model by 
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allowing only a single test of each variable. The restriction of polynomial width of quantum OBDDs 
leads to the number of qubits which is logarithmic in the size of the input. On the other hand, the gen- 
eralized lower bound on the width of quantum OBDDs [2 J leads to logarithmic number of qubits as a 
lower bound for almost all Boolean functions in the OBDD setting. 

For the model of quantum OBDDs we develop the fingerprinting technique introduced in 0. The 
basic ideas of this approach are due to Freivalds (e.g. see the book lfT3l ). These ideas were later success- 
fully applied in the quantum automata setting by Ambainis and Freivalds in 1998 [6] (later improved in 
0). Subsequently, the same technique was adapted for the quantum branching programs by Ablayev, 
Gainutdinova and Karpinski in 2001 [1], and was later generalized in 0. 

For our technique we use a variant of polynomial presentation for Boolean functions, which we 
call characteristic. The polynomial presentations of Boolean functions are widely used in theoretical 
computer science. For instance, an algebraic transformation of Boolean functions has been applied in 
iPTOl and for verification of Boolean functions. In the quantum setting polynomial representations 
were used for proving lower bounds on communication complexity in as well as for investigating 
query complexity in iTTSl . Our approach combines the ideas similar to the definition of characteristic 
polynomial from tfTOl . and to the notion of zero-error polynomial (see, e.g. lfl"8l ). 

Finally, we use the technique of polynomial projections to outline the limits of our fingerprinting 
method. This technique was intensively applied by Sauerhoff in the model of classical branching pro- 
grams (see, e.g., lfl5TD . In this paper we apply this approach for the quantum OBDD model by showing 
that the class of functions effectively computable via the quantum fingerprinting technique is closed 
under polynomial projections. 

2 Preliminaries 

We use the notation \i) for the vector from "K d , which has a 1 on the i-th position and elsewhere. An 
orthonormal basis |1),. . . ,\d) is usually referred to as the standard computational basis. In this paper we 
consider all quantum transformations and measurements with respect to this basis. 

Definition 1. A Quantum Branching Program Q over the Hilbert space "K d is defined as 

Q = {T,\y Q ), Accept), 

where T is a sequence of I instructions: Tj = [Xi-,Uj(0),Uj(l)J is determined by the variable Xj. tested 
on the step j, and Uj(0), Uj(l) are unitary transformations in l K d . 

Vectors |y/) 6 "K d are called states (state vectors) of Q, |y/o) £ "K d is the initial state of Q, and 
Accept C { 1 , 2, ... d} is the set of indices of accepting basis states. 

We define a computation of Q on an input G = G\ . . . G n G {0, 1}" as follows: 

1. A computation of Q starts from the initial state \ \jfo); 

2. The j-th instruction of Q reads the input symbol G{. (the value ofxij) and applies the transition 
matrix Uj = Uj(Oij) to the current state to obtain the state = Uj(Ojj) \ 

3. The final state is 
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4. After the l-th (last) step of quantum transformation Q measures its configuration = 
(0Ci , . . . , OLd) T , and the input O is accepted with probability 

Pr accept {p)= £ |a,-| 2 . 

iEAccept 



Note, that using the set Accept we can construct M accept - a projector on the accepting subspace 
"^■accept (i- e - a diagonal zero-one projection matrix, which determines the final projective measurement). 
Thus, the accepting probability can be re-written as 

Praccept{G) = (¥aMl ccept \M accept Ya) = \\Maecept |Va) Ill- 
Note also that this is a "measure-once" model analogous to the model of quantum finite automata 
in 021, in which the system evolves unitarily except for a single measurement at the end. We could also 
allow multiple measurements during the computation, by representing the state as a density matrix p, 
and by making the Uj superoperators, but we do not consider this here. 



Circuit representation. Quantum algorithms are usually given by using quantum circuit formalism 
||9j[T9l, because this approach is quite straightforward for describing such algorithms. 

We propose, that a QBP represents a classically-controlled quantum system. That is, a QBP can be 
viewed as a quantum circuit aided with an ability to read classical bits as control variables for unitary 
operations. 



=0= 



=0= 



IVo) 



r i0i> 
\<k) 



Ui(0) 



U 2 (l) 



V Z (0) 



=0= 



Ui(0) 



Here x,-, , . . . ,x, ; is the sequence of (not necessarily distinct) variables denoting classical control bits. 
Using the common notation single wires carry quantum information and double wires denote classical 
information and control. 



Complexity measures. The width of a QBP Q, denoted by width(Q), is the dimension d of the corre- 
sponding state space < K d , and the length of Q, denoted by length(Q), is the number / of instructions in 
the sequence T. There is one more commonly used complexity measure - the size of Q, which we define 
as size(Q) = width(Q) ■ length(Q). 

Note that for a QBP Q in the circuit setting another important complexity measure explicitly comes 
out - a number of quantum bits, denoted by qubits(Q), physically needed to implement a corresponding 
quantum system with classical control. From definition it follows that \ogwidth(Q) < qubits(Q). 
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Acceptance criteria. A QBP Q computes the Boolean function f with bounded error if there exists an 
£ G (0, 1/2) (called margin) such that for all inputs the probability of error is bounded by 1/2 — e. 

In particular, we say that a QBP Q computes the Boolean function f with one-sided error if there 
exists an £ G (0, 1) (called error) such that for all a G / _1 (1) the probability of Q accepting a is 1 and 
for all a G f~ l (0) the probability of Q erroneously accepting a is less than e. 

Read-once branching programs. Read-once BPs is a well-known restricted variant of branching pro- 
grams irm . 

Definition 2. We call a QBP Q a quantum OBDD (QOBDDj or read-once QBP if each variable x G 
{jci ,x n } occurs in the sequence T of transformations of Q at most once. 

For the rest of the paper we're only interested in QOBDDs, i.e. the length of all programs would be n 
(the number of input variables). Note that for OBDD model size(Q) = n • width(Q) and therefore we're 
mostly interested in the width of quantum OBDDs. 

Generalized Lower Bound. The following general lower bound on the width of QOBDDs was proven 



Theorem 1. Let f{x\ ,x n ) be a Boolean function computed by a quantum read-once branching pro- 
gram Q with bounded error for some margin £. Then 



where P is a deterministic OBDD of minimal width computing f{x\ > . • • >x n ). 

That is, the width of a quantum OBDD cannot be asymptotically less than the logarithm of the width 
of the minimal deterministic OBDD computing the same function. And since the deterministic width of 
many "natural" functions is exponential [ 17"], we obtain the linear lower bound for these functions. 

Let bits{P) be the number of bits (memory size) required to implement the minimal deterministic 
OBDD P for / and Q is an arbitrary quantum OBDD computing the same function. 

Then theorem [T]can be restated as the following corollary using the number of bits and qubits as the 
complexity measure. 

Corollary 1. 



3 Algorithms for QBPs Based on Fingerprinting 

Generally [13], fingerprinting - is a technique that allows to present objects (words over some finite 
alphabet) by their fingerprints, which are significantly smaller than the originals. It is used in randomized 
and quantum algorithms to test equality of some objects (binary strings) with one-sided error by simply 
comparing their fingerprints. 

In this paper we develop a variant of the fingerprinting technique adapted for quantum branching 
programs. At the heart of the method is the representation of Boolean functions by polynomials of 
special type, which we call characteristic. 



in H. 



width(e) > 



log width (P) 




qubits(Q) = Q.(logbits(P)). 
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3.1 Characteristic Polynomials for Quantum Fingerprinting 

The following definition is similar to the algebraic transformation of Boolean function from [ 10 1 and 0, 
though it is adapted for the fingerprinting technique. 

Definition 3. We call a polynomial g(x\ , . . . ,x n ) over the ring Z m a characteristic polynomial of a 
Boolean function f(x[ ,x n ) and denote it gf when for all o 6 {0, 1}" gf{o) = iff f(o) = 1. 

Lemma 1. For any Boolean function f there exists a characteristic polynomial gf over %2 n - 

Proof. One way to construct such characteristic polynomial g / is transforming a sum of products repre- 
sentation for -if. 

Let Ki V . . . V Ki be a sum of products for ->f and let Kj be a product of terms from K (negations ->Xj 
are replaced by 1 —xf). Then K\ + . . . + K is a characteristic polynomial over Z2» for / since it equals 
<J=^ all of Kj (and thus Ki) equal 0. This happens only when the negation of / equals 0. □ 

Generally, there are many polynomials for the same function, but we're interested in having a linear 
polynomial if it exists. 

For example, the function EQ„, which tests the equality of two ra-bit binary strings, has the following 
polynomial over : 

n n 

foC 1 -yd + (! - x dyd = L (•*« +?» ~ 2x '^) • 

(=1 !=1 

On the other hand, the same function can be represented by the polynomial 



v'-i 

1=1 1=1 



Some functions don't have a linear characteristic polynomial at all (e.g. a disjunction of n variables 
/ = jci VX2 V . . . Vi„), while their negations have a linear characteristic polynomial (e.g. g^f = Yli=\ x i 
over Z„ + i). 

We use this presentation of Boolean functions for our fingerprinting technique. 
3.2 Fingerprinting technique 

For a Boolean function / we choose an error rate £ > and pick a characteristic polynomial g over 
the ring Z OT . Then for arbitrary binary string a = 0\...o n we create its fingerprint \h a ) composing 
? = 2 riog((2/ e )in2m)l s i ng i e qubit fingerprints \h' a ): 

\h a ) = cos 0) +sin —- 1) 

into entangled state \h a ) of log? + 1 qubits: 

\ho) = ^tW°)- 
v t i=l 

Here the transformations of the last qubit in t different subspaces "simulate" the transformations of all of 
the Wa) = 1, • • • ,0- That is, the last qubit is in parallel rotated by t different angles about the y axis of 
the Bloch sphere. 

The chosen parameters G {1, . . . ,m— 1} for i G {1, . . . ,t} are "good" following the notion of [6]. 
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Definition 4. A set of parameters K = {k\ , . . . , k t } is called "good" for some integer b ^ mod m if 

1 /' 2nkib\ 

j2 1 L cos — - 1 f - 



m 



The left side of inequality is the squared amplitude of the basis state |0)® logf |0) if b = g(o) and the 
operator //® logf <g)/ has been applied to the fingerprint \h a ). Informally, that kind of set guarantees, that 
the probability of error will be bounded by a constant below 1 . 

The following lemma proves the existence of a "good" set and generalizes the proof of the corre- 
sponding statement from Q. 

Lemma 2. ffi There is a set K with \K\ = t = 2 r io s(( 2 / e ) ln2 ' OT )l w hich is "good" for all integer b / 
mod m. 

We use this result for our fingerprinting technique choosing the set K = {k\ ,...,k t } which is "good" 
for all b = g(o) 7^ 0. That is, it allows to distinguish those inputs whose image is modulo m from the 
others. 



3.3 Boolean Functions Computable via Fingerprinting Method 

Let f{x\ , . . . ,x n ) be a Boolean function and g(x\ ,x n ) be its characteristic polynomial. The following 
theorem holds. 

Theorem 2. Let e G (0, 1). If g is a linear polynomial over 7L m , i.e. g = c\X\ + ... c n x n + cq, then f can 
be computed with one-sided error e by a quantum OBDD of width O 

Proof. Here is the algorithm in the circuit notation: 

X\ 



t - 



H 



|01ogf) -\H_ 
\§ tea-get) 



-O 



-o 



-6 



RlA Rf.n ^1,0 Rt.O 



H -M- 



H 



t t t 

IVo) IVi) |V2> 



t 



t t 



Initially qubits ® |0 2 ) <8>- • -® |<frogj)® garget) are in the state |^ ) = |0)® log ' |0). For i G {1, . . . ,t}, 
j G {0, ...,«} we define rotations as 



AnkjCj 



m 



where cj are the coefficients of the linear polynomial for / and the set of parameters K = {k\ ,k t } is 
"good" according to the Definition g] with t = 2riog((2/e)ln2-m)l_ 
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Let a = (7i . . . o n 6 {0, l}' 1 be an input string. 

The first layer of Hadamard operators transforms the state | \ffo) into 

|Vl) = 4rI>'>|0>- 

V' j=l 

Next, upon input symbol identity transformation / is applied. But if the value of xj is 1, then 
the state of the last qubit is transformed by the operator Rjj, rotating it by the angle proportional to Cj. 
Moreover, the rotation is done in each of t subspaces with the corresponding amplitude l/y/i. Such 
a parallelism is implemented by the controlled operators Cj(Rij), which transform the states \i) |-) into 
\i)Rij |-), and leave others unchanged. For instance, having read the input symbol x\ = 1, the system 
would evolve into state 

|V2> = 4 ? EQ(/? i) i)|i)|0) = 4.E|i)/? i)1 |0) 
v i=i v i=i 

= ^Ll')(cos^|0) + sin^a|l)) ' 
Thus, after having read the input o the amplitudes would "collect" the sum £" =1 cyO/ 

IW = jg E |«> (cos 'V J J |0)+sin 'V I 1 )) ■ 

i=i v 7 

At the next step we perform the rotations by the angle 47r ^ c ° about the y axis of the Bloch sphere for 
each i € {1, . . . ,*}. Therefore, the state of the system would be 

IV4> = ^IJ/)(cos^|0) + sin^l|l)). 
Applying // 8l *®/we obtain the state 

|V5> = (|lcos^l)|0)® logf |0) + 

+y|0)® to 8 r |l>+EI'")(a,-|0) + ft|l)) ) 

where 7, a,, and j3, are some unimportant amplitudes. 

The input a is accepted if the measurement outcome is |0)® log ' |0). Clearly, the accepting probability 

is 

p („\ 1 / v g^g(g) V 

Praccept(o) = 72 I L COS ^ J ' 

If /(a) = 1 then g(a) = and the program accepts a with probability 1. Otherwise, the choice of 
the set K = {ki , . . . , k t } guarantees that 

. . I f 4, 2nkig{o)\ 2 

Pr a ccept{0) = 72 ( L COS J < £ ' 

Thus, / can be computed by a quantum OBDD Q, with qubits(Q) = log? + 1 = 0(log (^j^j )■ The 
width of the program is 2<i Mt -<Q) = O (^f^j . □ 
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The following functions (for definitions see, e.g., E)) have the aforementioned linear polynomials 
and thus are effectively computed via the fingerprinting technique. 



MOD m 


L Xi 

1=1 




0(loj 


;logra) 


MUU m 


n 
1=1 


/L, m 


C(iO£ 


; log m) 


EQn 


;=i i=i 




0(lo£ 


;«) 


Palindrome n 


L»/2J 

1=1 i=[»/2l 




0(loj 


I") 


PERM n 


1 1 xijdn+iy-i+in+iy+j- 1 ) 

1=1 7=1 
2n 

1=1 




0(loj 













4 Generalized Approach 

The fingerprinting technique described in the previous section allows us to test a single property of the 
input encoded by a characteristic polynomial. Using the same ideas we can test the conjunction of several 
conditions encoded by a group of characteristic polynomials which we call a characteristic of a function. 

Definition 5. We call a set Xf of polynomials over Z m a characteristic of a Boolean function f if for all 
polynomials g G Xf an d a ^ ® £ {0, 1}" it holds that g(o) = iff G G 

We say that a characteristic is linear if all of its polynomials are linear. 

From Lemma [Hit follows that for each Boolean function there is always a characteristic consisting 
of a single characteristic polynomial. 

Now we can generalize the Fingerprinting technique from section [3721 

Generalized Fingerprinting technique For a Boolean function / we choose an error rate e > and 
pick a characteristic Xf = {gi> • • • ,8l} over ^m- Then for arbitrary binary string o = Ci . . . o n we create 
its fingerprint \h a ) composing t -l (t = 2r io s(( 2 / e ) ln2m )l) single qubit fingerprints |/i' ff (./)): 

KU)) = 

K) = 

Theorem 3. If Xf ' s a linear characteristic then f can be computed by a quantum OBDD of width 
0(2l*/l log m). 

Proof. The proof of this result somewhat generalizes the proof of Theorem |2]to the case of several target 



cos n ' 0) +sin n ' 1) 

^ L\i) |^(1)> |^(2)>- |*U0) 
i=i 
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qubits. Here is the circuit: 



|<M -\H\ $ 

102) 



|01ogf> 
\$ target ) 

k 2 v 



fii> 



-0 



1.1 



-o 



ft*) 



-Q 



-O 



Ml) 



ft 





ft?, 















ft 



1.0 



- ft^l | 1 ^1,0 I 



R 



ft 



1.0 





1 

.0 




R 


2 
.0 






R t,0 



1. Upon the input a = 0\ . . . o n we create the fingerprint \h a ). 

2. We measure \h a ) in the standard computational basis and accept the input if the outcome of the 
last / qubits is the all-zero state. Thus, the probability of accepting a is 



Pr, 



1 



accept (<*) = tI^ cos ' 

1 i=\ 



7ik igl (a) 



•cos 



2 nkigiia) 



m 



m 



If f(o) = 1 then all of gi(o) = and we will always accept. 

If f(o) = then there is at least one such j that gj(o) ^ and the choice of the "good" set K 
guarantees that the probability of the erroneously accepting is bounded by 



Pr, 



accept 



\ t cos 2 



2 ickigija) _ _ _ ? xkjgi(o) 
m 

i=i 

2 xkgjjo) _ J_ y 



< 7 I cos 

i=i 



in 



iX,i(l + cos 



2 + 2? I COS 



2nkigj{a) 



2 ' 



The number of qubits used by this QBP Q is qubits{Q) = 0(loglogm + \%f\). Therefore, the width 
of the program is 2" ubits( -® = 0(2)Xf\ logm). 

□ 

Note that though this upper bound is exponential our approach can be effectively used when the size 
of a characteristic is O(loglogm) and m = 2 n ° W . That is, in this case we will stay within the polynomial 
width. 

The Theorem [3] has two immediate consequences which might be useful for proving upper bounds. 
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Corollary 2. If a Boolean function f = f\ A/2 A . . . Af s is a conjunction ofs = 0(log n) Boolean functions 
fi, each having a linear characteristic polynomial over 7Li n > then f can be computed by an 0(\ogn)-qubit 
quantum OBDD. 

Corollary 3. If a Boolean function f = f\ V/2 V . . . Vf s is a disjunction ofs = 0(\ogn) Boolean functions 
fi and the negation of each has a linear characteristic polynomial over Z2*, then f can be computed by 
an 0(logn)-qubit quantum OBDD. 

The last corollary uses the fact that -if = -1/1 A -1/2 A ... A -<f s . 

The generalized approach can be used to construct an effective quantum OBDD for the Boolean 
variant of the Hidden Subgroup Problem H. 

5 Using reductions for fingerprinting 

Reduction is a well-known concept in complexity theory. The most investigated reduction for the model 
of branching programs is the technique of polynomial projections |fT7|. 

Definition 6. The sequence (/„) of Boolean functions is a polynomial projection of(h n ), (f n ) <p ro j {hn), 
if fn = hp(n)(yi:' ■ ■ ■>yp{n))f or some polynomial p and yj £ {0,1, ,x\, . . . ,x n ,x n }. The number of j such 
that yj E {jtj,x,} is called the multiplicity ofxj. 

Note that this type of reduction keep the linearity of the corresponding characteristic polynomials 
and thus the effective computability via the fingerprinting method. 

Lemma 3. If (/„) <^> ro j (h n ) and each h n can be represented by a linear characteristic, then each f n 
also has a linear characteristic. 

Proof. The variable substitutions corresponding to projection are y, = 0, y, = 1, y\ = Xj t ,yi = 1 —Xj r Ob- 
viously, if a characteristic is linear then substitutions of these types to polynomials g = c\y\ + ... c n y n + co 
keep them linear. □ 

On the other hand, the a special case of polynomial projections can be used to prove lower bounds in 
the OBDD model. 

Definition 7. A projection is called a read-once projection, (/„) <rop (h n )> if the multiplicity of each 
variable is bounded by one. 

The relation <rop is reflexive and transitive, moreover if (/„) <rop (h n ) and the OBDD size of h n 
is bounded by the polynomial q(n) then the OBDD size of /„ is bounded by the polynomial q(p(n)). 
The last property was proved for the classical OBDD model (see, e.g. ifTTl "). but it can be proved for the 
quantum setting in an analogous way. 

Thus the read-once projections may be used for proving lower bounds of the subclasses of Boolean 
functions whose projections are exponentially hard for the OBDD model. These are, for example, Set- 
Disjointness and Neighbored Ones Boolean functions which were proved to be exponentially hard in 

na. 

In particular these subclasses of "hard" functions cannot be effectively computed via quantum fin- 
gerprinting method in the model of read-once quantum branching programs. 

Overall, the technique of polynomial projections outlines the limits of our fingerprinting method in 
the quantum OBDD model with upper bounds propagated by the general polynomial projections and 
lower bounds based on individual "hard" functions and read-once projections. 
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